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Graphene and its quasi-one-dimensional counterpart, graphene nanoribbons, present an ideal platform for 
tweaking their unique electronic, magnetic and mechanical properties by various means for potential 
next-generation device applications. However, such tweaking requires knowledge of the electron-electron 
interactions that play a crucial role in these confined geometries. Here, we have investigated the magnetic 
and conducting properties of zigzag edge graphene nanoribbons (ZGNRs) using the many-body 
configuration interaction (CI) method on the basis of the Hubbard Hamiltonian. For the half-filled case, the 
many-body ground state shows a ferromagnetic spin-spin correlation along the zigzag edge, which supports 
the picture obtained from one-electron theory. However, hole doping reduces the spin and charge excitation 
gap, making the ground state conducting and magnetic. We also provide a two-state model that explains the 
low-lying charge and spin excitation spectrum of ZGNRs. An experimental setup to confirm the 
hole-mediated conducting and magnetic states is discussed. 
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Carbon nanomaterials have gained sustained interest in recent times due to their fascinating electronic 
properties, which arise from electronic confinement in a reduced-length scale. The successful isolation of 
graphene by mechanical exfoliation' has provided further impetus and has allowed for the first time the 
understanding of various properties in a truly two-dimensional contexf The ultrahigh charge carrier mobility, 
transparency and mechanically flexible properties of graphene provide an excellent platform for futuristic device 
applications'. 

Another intriguing aspect of graphene is the strong nanoscale and edge effects on its electronic and magnetic 
properties. The quasi-one-dimensional ZGNRs have drawn particular attention because of their peculiar edge 
ferrimagnetism that arises from the edge-localised states near the zigzag edges" '°. Intensive studies ranging from 
mean-field to density functional theory have been performed to understand the carbon-based magnetism in 
ZGNRs" A recent experiment using scanning tunnelling spectroscopy (STS) reports the splitting of the density 
of states due to the edge magnetism for chiral graphene nanoribbons'^. 

These theoretical and experimental studies have motivated the fabrication of spintronics devices based on 
graphene nanoribbons. Although its edge-state magnetism seems to contradict the Mermin- Wagner theorem, 
which rules out the possibility of long-range ordering in quasi-one-dimensional systems"", a recent theoretical 
report using quantum Monte Carlo supports the existence of such a long spin-spin correlation length along the 
zigzag edge and justifies the picture obtained from the one-electron theories". Moreover, many previous studies 
suggest various ways to tweak the electronic properties of ZGNRs to achieve magneto-transport by means of 
doping, chemical modifications or external perturbations'""^^. 

From the viewpoint of device applications, it is necessary to clarify the interplay between the edge magnetism 
and the hole-doping effect, as the electron density in the graphene system can be easily tuned using the back gate 
electrode^^'"^^. However, the hole-doping effect on edge magnetism has not been studied with the appropriate 
inclusion of electron-electron interactions. 

In this study, we theoretically investigate the magnetic and conducting properties of ZGNRs and their 
response to doping. The development of gapless charge and spin excitations with hole doping, i.e., hole- 
mediated metallic ferromagnetism, is also discussed for the first time in terms of its potential applications. 
Our theoretical analysis is based on a large-scale numerical simulation using the many-body configuration 
interaction (CI) method with the complete active space (CAS-CI) approximation, which correctly includes 
quantum fluctuations. We show that the microscopic origin of magnetism in ZGNRs can be well understood 
on the basis of a generic two-state model. 
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Results 

We model the systems within a Hubbard Hamihonian with onsite 
electron -electron interactions. The Hamiltonian reads as 



H=-t 



where a]^ creates (annihilates) one electron with spin a on site / 
and tiju is the number operator. "U" denotes the onsite Coulomb 
interaction, and the nearest neighbour hopping integral "t" has been 
set to unity. We restrict our calculations within the weak coupling 
limit and consider the U/t ratio as one. This assumption ensures the 
absence of polarisation of spins for the bulk sites in the graphene. 
Note that the above Hamiltonian conserves the electron-hole sym- 
metry, and therefore, hole doping and electron doping are expected 
to show identical behaviour. 

Figure la shows a schematic representation of the ZGNRs. The 
lattice constant has been considered to be unity. The width of the 
ZGNRs is defined by the number of zigzag lines (Nj.) along the 
transverse ribbon direction. The edges are assumed to be passivated 
by hydrogen atoms, and therefore, the electronic states near the 
Fermi energy remain unaltered. Within the non-interacting case, 
the top of the valence band and the bottom of the conduction band 
are partially flat between k = n (the Brillouin zone (BZ) boundary) 
and k = 2Tt/3 (the Dirac point) at the Fermi energy. As shown in 
Figure lb, this flat region increases with increasing ribbon width. 

We calculate the many-body charge gap (Ac) and spin gap ( A^) and 
their response to the hole doping as foUows^*"'^': 

A,{k) = [Eo{K + l,So) +Eo{N,-l,So)-2Eo{N„So)]{k) 



As{k) = [£o(N„So-h 1) -£o(N„So)](/:) 

Here, Eo(N^So) represents the lowest energy in the subspace of 
( = N— Hf,) electrons, where n^ denotes the number of holes doped 
in the half-filled system (N) with the lowest SJ" (the z-component of 
total spin) = So ( = 0 or 1/2 for even or odd N^, respectively) at each k- 
point. In this paper, to probe the different spin sectors, we vary the 
Sz'"'. The number of electrons and SJ" both commute with the 
Hubbard Hamiltonian and therefore can be used as good quantum 
numbers. Note that the A^ and A^ are calculated at discrete fc-points, 
and the scattering processes connecting the different fc-points are 
neglected. This approximation occasionally fails to capture the essen- 
tial physics of spin excitations in the case of narrow ribbons, which 
will be discussed in the spin gap section. However, the low energy 
charge and spin excitation properties for wider ribbons can be sat- 
isfactorily captured by this approximation, as is evident from our 
calculations for finite size systems in the latter part of the article. 

Two-state model. In a weak coupling regime, the charge and spin 
excitation properties of the ZGNRs are governed by two bands 
(shown in Fig. lb). Thus, the whole problem can be effectively 
mapped into a two-state model for the half-filled case. At each k- 
point, the two-state model is defined by four possible basis states, as 
shown in Figure Ic. In the (7=0 limit, the system is stabilised only by 
kinetic energy, and the ground state is defined only by ^j. 

Before showing the detail of the two-state model, we shall try to 
avoid confusion between molecular orbital space and real space 
descriptions, as the two-state model is constructed in molecular 
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Figure 1 | Geometry and band structure of ZGNRs with two-state model, (a) The structure of ZGNR, with being the width. The unit cell (dashed 
rectangle) consists of same number of A (filled circles) and B (open circles) sublattice points, (b) Tight binding energy spectrum of ZGNR in momentum 
(k) space. The top of the valence band and the bottom of the conduction band of ZGNRs with JVj. = 4 and 20 are depicted as the solid lines. These two 
central subbands construct the two-state model. The shaded region designates the bulk graphene spectrum near the Fermi energy (£ = 0) with linear 
dispersion around Dirac point (k = 2n/3). (c) The four basis states ( ijdj, (^21 '(''3 and (^4) to construct the two-state model for half-filled system with 5;.'°' = 0. 
(d) W(N^So) is the wave function for the ground state for half- filled system. *P{N^ + l,So) and W(N^ - hSo) are the wave functions for the system with one 
more and one less electron than half-fiUing, respectively. ¥'(N„So + 1) represents the ground state of half-fiUed system with SJ"' = 1. The wave function 
coefficients are represented in terms of c. 
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Figure 2 | Charge gap, spin gap and spin density of half-filled ZGNRs. (a) The charge gap (Ac) behavior over the BZ for ZGNRs with Nj. = 4 and 20 (solid 
lines) and for bulk graphene (dashed line) at half-fiUing. (b) The minimum of charge gap, (min.) (circles) near the Dirac point as a function of and 
fitted (solid line) with the equation: Ac= 1.736/ N^. The zl^ at A: = n for varying is shown by dashed line with squares, (c) The spin gap (zlj behavior 
over the BZ for ZGNRs with = 4 and 20 (solid lines) and for bulk graphene (dashed line) at half-filling, (d) and (e) shows the spin density profile of 
half-filled ZGNR with Nj, = 20 for SJ"' = So and So + 1, respectively as a function ofwavenumber (k) along the width (designatedby site index). The color 
bar shows the absolute scale of spin density. 



orbital space rather than in real space. Let us first discuss the simplest 
example of a system with two atomic sites. As we know, in real space, 
the introduction of U begins localising one electron at each atomic 
site. Therefore, the ground state for half-filling and 5^'°' = 0 with 

non-zero L/can be represented as a singlet state, —= 1 1 J. ~ J.t)> in real 

V 2 

space. In molecular orbital space, this ground state can be described 

as 101 — ^4), which corresponds to 'P(N^,So) with equal Cj and C2, 
V2 

as shown in Figure Id within the two-state model. Note that (^2 and 
03 contribute to the excited states but not to the ground state. 

The charge excitation process can be understood based on the two- 
state model as follows. The effect of 1/ becomes more pronounced over 
the flat band region (27i/3 < k s ti for wide ZGNRs), as the on-site 
Coulomb interaction, U, is effectively evaluated as U/W, where W is 
the bandwidth^**. Because W is very small over the flat band region for 
any U > 0, the effective interaction becomes very high. Thus, the 
groimd state, 'P(N^So), is described as a Unear combination of 
and 04 (see Fig. Id), which has been confirmed by analysing the 
coefficients of the wavefimctions obtained from the CI calculations. 
However, as we move away from the flat band region (|kj s 2n/3), the 
decrease in U/W increases the contribution of 0j to the ground state. 

In the case of one fewer electron than half-fiUing, the ground state 
is singly occupied in the lower flat band and is expressed as ^(N^ - 
l,So), as shown in Figure Id. Due to the single occupancy, this state is 
independent of the parameter U. However, the state with one more 
electron than half-filling cannot avoid the pairing of electrons, and 
hence the ground state, 'P(Ni. + l,So), costs 17 amount of energy. 
Therefore, the difference between EgiN^ +l,So) and Eo(Ni,— l,So) 
always remains JJ"'. Thus, the charge gap, A^, is solely determined 
by the relative position of EgiN^So) with respect to the states with one 
extra and one fewer electron. This behaviour will be confirmed in the 
CAS -CI calculation later. 

The two-state model provides a microscopic understanding of the 
spin gap behaviour within the many-body framework. Because in the 



spin excited state, W(Ng,So +1) (shown in Fig. Id), the valence and 
conduction bands are occupied by single electrons with the same 
spin, its energy Eo(Ne,So +1) becomes independent of U. Thus, A^ 
is determined only by Eo(Ni,, So)- The most intriguing behaviour of the 
zero spin gap is observed when 'P(Ng,So) has equal contributions 
from (f>i and 0^. The top of the valence band is doubly occupied in 
the 0j configuration, whereas in the 0^ configuration, the bottom of 
the conduction band is doubly occupied (see Fig. Id). When each of 
them contributes 50% to construct the ground state, it effectively 
becomes similar to the state with single occupancy of the top of the 
valence band and the bottom of the conduction band. Therefore, the 
energy of the ground state, EgiN^Sg), becomes equal to the energy of 
the spin excited state, Eo(N^,So +1), closing the spin gap. This cri- 
terion actually can be satisfied in the flat band region, where U/W is 
large enough. Because of the vanishing W over this region, any U> 0 
is enough to make the ground state a magnetic state. This magnetic 
ground state, arising from degenerate flat bands near the Fermi 
energy, is analogous to Mielke's flat band ferromagnetism^''"*. 

Charge gap. We characterise the conducting properties of ZGNRs by 
evaluating the charge gap, A„ which is the energy difference between 
charging and discharging an electron to the system. The decrease in 
the charge gap indicates the enhancement of conduction behaviour. 
Our theoretical analysis supports the pictures obtained from one- 
electron theory regarding charge excitation. 

Figure 2a shows the charge gap A^ for half- filled ZGNRs over the 
BZ using CAS-CI. The observations within the two-state model 
remain unchanged even for higher numbers of bands in the CAS- 
CI calculations, where 4900 numerical basis states were considered. 
The minimum in A^ decreases and shifts towards the Dirac point 
with increasing width to resemble the conical spectrum of bulk gra- 
phene. The minimum charge gap. Admin.), varies inversely with the 
width of the ZGNRs (see Fig. 2b) and approaches zero at infinite 
width, i.e., bulk graphene. All of these observations are qualitatively 
in agreement with the one-electron approaches such as mean-field or 
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local density approximation (LDA)'^'^' and with the analytic 
results'". Interestingly, in this weak coupling regime, the at BZ 
boundary (fc = 7t) is exactly equal to the U/t ratio and remains 
unaltered with increasing width (see Fig. 2b). This equality originates 
from the fact that at A: = 7t, the ground state energy has a vanishing 
contribution from the kinetic energy term of the Hubbard 
HamUtonian and is governed only by the interaction parameter. In 
the weak coupling regime, the half-filled system and the system with 
one fewer electron both become the ground state of the interacting 
part of the Hamiltonian, and thus, the energy of removing one elec- 
tron from the half-filled state amounts to zero at A: = 7i, resulting in a 
charge gap of U/t. 

Spin gap. To investigate the behaviour of spins in ZGNRs, we 
calculate the spin gap, A^, which is the amount of energy required 
to flip one spin in the system. With decreasing A^, the system 
becomes more subject to magnetic excitations. A zero spin gap 
refers to the degenerate magnetic ground state. 

Figure 2c shows the A^ for half-filled ZGNRs of different widths 
over the BZ, obtained using the CAS-CI with 4900 basis states. With 
increasing width, the ZGNRs show gapless spin excitation over a 
larger region within k = n and the Dirac point. This behaviour is 
due to the increasing flat band region with increasing width, as we 
have described within the two-state model. Note that the appearance 
of a zero spin gap is unlikely for very narrow ribbons with only two or 
three zigzag chains, as evident from previous studies on spin lad- 
^gj.gii,32 ji^g difference in spin excitation behaviours for odd and 
even numbers of legs of spin ladders is not observed in our calcula- 
tions'*'^'''. This discrepancy can be attributed to the absence of back- 
ward scattering and umklapp scattering in our calculations. 
Nevertheless, we wUl present the systematic investigation of finite 
size systems, including these scattering processes, later in this article 
to validate our observations. Moreover, it has been observed that the 
spin gap progressively decreases with increasing width of the ladder. 
Therefore, we emphasise that our observations are robust in the case 
of wider ZGNRs. 

The antiparallel alignment of spins within the opposite edges of 
ZGNRs'*''''^''''"''"' remains unaltered within the many-body level of 
theory, due to the bipartite nature of the graphene lattice. As seen 
in Figure 2d, atk = n, the spins are completely localised in antipar- 
allel fashion on the opposite edges for half-filled ZGNRs with SJ"' = 
So- However, a gapless spin excitation over the flat band region indi- 
cates that the state with Sz'°' = So + 1, which shows parallel spin 
orientation between opposite edges (see Fig. 2e), is degenerate to the 
antiferromagnetic ground state, as suggested by Lieb's theorem'^. 
This observation suggests that once the complete ferromagnetic state 
is achieved by applying an external magnetic field, it persists as a 
stable state even after switching off the magnetic field in the case of 
undoped ZGNRs. 

Hole doping. Now, let us investigate the magnetic and conducting 
properties of ZGNRs in the presence of holes. Our calculation 
indicates the reduction of charge and spin gaps upon hole doping, 
i.e., hole-induced metallic magnetism. Figure 3a shows the charge 
gap A^. variation with hole doping over the BZ for = 20 ZGNR. 
The charge gap decreases drastically as the system moves away from 
the half-filled state, indicating a hole-induced transition towards the 
metallic phase. Note that the A^ for even and odd numbers of holes 
merge to two different values near the Dirac point, due to the 
systematic removal of electrons from the doubly occupied valence 
bands. However, this odd-even effect becomes indistinguishable for 
ZGNRs with larger widths when the charge gap for the hole-doped 
ZGNRs becomes negligibly small. Therefore, we observe a hole- 
mediated opening of conduction channels through the ZGNRs. 

Figure 3b shows the spin gap A^ for different hole dopings for = 
20 ZGNR. For single-hole doping, the spin needs to be excited from 
the band below the top of the valence band to the bottom of the 
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Figure 3 | Charge gap and spin gap of hole-doped ZGNRs. (a) The charge 
gap {Ac) and (b) spin gap (A^) of ZGNR with = 20 over k = nl2 n for 
different hole (n^) doping. 

conduction band. Therefore, the spin gap in this case is always higher 
than in the half-filled case, where the spin excitation comes from the 
top of the valence band. However, with increasing width, the closer 
spacing of the two top valence bands makes the spin gap for both 
cases almost equal over the k region from the BZ centre (fc=0) to the 
Dirac point. However, from the Dirac point to the BZ boundary, the 
zlj increases due to the appearance of edge states. The A^ for systems 
up to three-hole doping also shows a monotonic increase over this 
region, due to spin excitation to the top of the valence band, i.e., the 
edge state. The spin gap for these systems with up to three holes, 
however, decreases over the k region from the BZ centre to the Dirac 
point. Higher hole doping results in a gradual closing of the spin gap 
over the full BZ, and consequently, the system turns out to he magnetic. 
With increasing width, the A^ decreases without altering the qualitat- 
ive behaviour. This decrease in the spin gap indicates hole-induced 
magnetism in ZGNRs. The observed odd-even effect in A^ also 
becomes indistinguishable with increasing ZGNR width. 

This appearance of magnetism drives us to explore the spin struc- 
ture in ZGNRs with hole doping. Figure 4 shows the momentum- 
dependent spin density profile over the full BZ for an = 20 ZGNR 
for the spin-excited state, i.e., SJ" = So + 1. With increasing hole 
doping, the spin-excited state of ZGNRs gradually stabilises and 
eventually turns out to be the ground state. The increase in width 
makes this process easier. As seen in the figure, the spin density 
mainly localises on the ZGNR edges at the BZ boundary for lower 
hole doping and gradually penetrates to the bulk sites as we move 
away from k = n. With increased hole doping, the absolute spin 
density on the edge atoms over the flat band region gradually 
decreases because of the systematic disappearance of electrons from 
the edge states. However, most importantly, all the remaining elec- 
trons align with the same spin orientation throughout the whole 
width, giving rise to a net magnetic moment for the hole-doped 
ZGNRs. 

For further verification, in Figure 5, we plot the spin-spin correla- 
tion ((Sf S/)) between two opposite-edge atoms {i and;') for a = 20 
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Figure 4 | Spin density profile of hole-doped ZGNRs. Spin density profile of hole doped ZGNR with N^= 20 and SJ"' = So + 1 as a function of 
wavenumber (k) along the width (designated by site index). The hole doping increases from left to right with (a) nh = 1, (b) nj, = 2, (c) nj, = 3 and (d) nj, = 
4. The color bar shows the absolute scale of spin density. 



ZGNR with SJ"' = Sg over the full BZ for different dopings. At fc = 7t, 
the half-flUed ZGNR shows complete antiferromagnetic correlation 
between the two opposite edges, which decreases monotonically 
towards the Dirac point. We observe that the correlation within 
the bulk sites increases as k-^lnli from its zero value at A: = Tt and 
that in the vicinity of Dirac point, all the correlations within the bulk 
and edge sites become comparable. 

The antiferromagnetic correlation between edges is largely 
destroyed by hole doping. Upon hole doping, electrons first vacate 
from the top of the valence band, i.e., the edge state, and hence, the 
spin density on the edge atoms vanishes gradually. In the case of an 
even number of holes, the spins disappear symmetrically from both 
the edges, but a system with an odd number of holes does not main- 
tain a symmetrical spin distribution. Therefore, the spin-spin cor- 
relation shows an odd-even effect (see Fig. 5), though it is negligibly 
small due to quantum fluctuations. We observe that with hole dop- 
ing, the antiferromagnetic correlation within two different sublattice 
points throughout the whole lattice also decreases over the full BZ. 
From the behaviour of the spin density profile and the spin-spin 
correlation, we can infer a hole-induced net magnetisation in 
ZGNRs. 

Thus far, we have neglected the effect of scattering events connect- 
ing two different fc-states to emphasise the simple and unified picture 
for the spin and charge excitation mechanism and their direct cor- 
respondence to the two-state model for ZGNRs. However, to invest- 
igate the effect of such scattering events, we perform the CAS-CI 
calculations for non-periodic finite size systems with long, smooth 
zigzag edges, which inherently include scattering processes arising 
from electronic correlations. 




Figure 5 | Spin-spin correlation of ZGNRs. Spin-spin correlation 
between two opposite edge atoms of ZGNR with = 20 over the BZ for 
different hole (nh) doping. The inset shows zoomed plot near (SfSf) = 0. 



Figures 6a and 6b present the charge-gap and spin-gap behaviours, 
respectively, for such finite size systems as a function of inverse 
length (L), i.e., the number of unit cells for a fixed width, = 6. 
Both zfj and decrease upon hole doping, indicating conducting 
and magnetic behaviours, respectively, as we have already seen for 
periodic ZGNRs (see Fig. 3a and 3b). Interestingly, the very low value 
of z<j even in the case of the half-fdled system (n^ = 0) suggests an 
easily accessible magnetic state, as we have observed over the flat 
band region of periodic ZGNRs (see Fig. 2c). 

We also investigated the spin-spin correlation behaviour for finite 
size systems with different hole dopings. Figure 6c shows the spin- 
spin correlation calculated from one edge atom (the "0" -th site in 
the finite graphene structure, shown in the middle panel) with all 
other edge atoms along the same edge (top panel) and with all other 
edge atoms along the opposite edge (bottom panel). The ferromag- 
netic correlation (positive (S/S/)) along the same edge decreases 
slowly with increasing hole doping because of the gradual disappear- 
ance of electrons from the edges, i.e., from edge states. The antifer- 
romagnetic correlation (negative {SgSf}) within the opposite edges 
shows a sharp decrease with hole doping, as observed in periodic 
ZGNRs (Fig. 5). All of these observations collectively lead us to the 
conclusion that the scattering events connecting different A:-points in 
momentum space do not alter the qualitative outcome of our study 
for periodic ZGNRs. 

Discussion 

The observation of the magnetic ground state of ZGNRs, along with 
the decreased charge gap, suggests that hole-doped ZGNRs can act as 
magnetic conductors with an abundance of one type of spin in the 
charge carriers. This behaviour can be exploited to fabricate spin 
filters by controlling the background dielectric, as shown schematic- 
ally in Figure 7. The hole doping can also be accomplished using 
oxidising agents that attract the electrons from the ZGNRs without 
destroying their backbone geometry. The originally antiferromag- 
netic and semiconducting ZGNRs (see Fig. 7a) can be turned into 
conducting spin filters by the application of gate bias (see Fig. 7b). 
However, although the device setup is rather conventionaP^'^'*"''^, 
experimentally achieving magneto-transport requires the graphene 
channel to have smooth zigzag edges and the other perturbations 
such as substrate effect and electrode coupling to be minimised. In 
such an experimental situation, spin-dependent transport might be 
detected at low temperatures. From a theoretical viewpoint, the 
knowledge of scattering process and velocity of the charge carriers 
is also necessary for a complete understanding of spin-filtering beha- 
viour. However, these issues are beyond the scope of the present 
study. Nevertheless, we look forward to further experimental and 
theoretical studies to realise this proposed experimental device setup. 
We are also pursuing transport calculations at the many-body level of 
theory in our follow-up studies. 
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Figure 6 | Charge gap, spingapandspin-spincorrelationof finite size graphene and their response to hole doping, (a) The charge gap (/(J and (b) spin 
gap (As) for finite graphene system with long zigzag edge and fixed — 6 as a function of inverse length (L), i.e., number of unit cells for different hole 
(n^i) doping, (c) The spin-spin correlation {(Sg'Sf)), calculated from the "0" -th atom, located at the middle of the edge of the graphene structure with 
— 6 (middle panel) with all the edge atoms (j) along the same edge (top panel) and with all the edge atoms (j) along the opposite edge (bottom panel) for 
different hole (Uh) doping. 



To summarise, we have investigated the electronic properties of 
ZGNRs and their response to hole doping at the many-body CI level 
of calculations, taking into account the quantum fluctuations. The 
charge and spin excitation processes in ZGNRs have been explained 
based on an effective two-state model, which considers the top of the 
valence band and the bottom of the conduction band. Our study 
suggests that hole doping not only causes the ZGNRs to become 
conducting but also reduces the spin gap to make the ground state 
magnetic. This magnetism can be enhanced further by creating nano 
pores with zigzag edges in the bulk of the ZGNRs. These observations 
indicate the possibility of fabricating a spin-filter device based on 
ZGNRs by tuning the doping level. The fundamental understanding 
of the charge and spin excitation behaviours will guide further 
research in graphene and related materials and will provide tremend- 
ous impetus to the fascinating research area of carbon-based mag- 
netism. 

Methods 

Because nanoscale materials contain a large number of atomic sites, it is very difficult 
to treat the electron-electron interactions properly. Although the mean field approach 
or the local density approximation (LDA) allows for a numerical solution of the 



Hubbard HamUtonian even for sufficiently large widths, the one-electron description 
of electronic correlations suffers severely from the lack of representation of the 
quantum fluctuations. However, the many-body methods such as the density matrix 
renormalisation group and exact diagonahsation are not applicable for larger system 
sizes because of dimensionality limitations^^. 

In this paper, we employ the many-body configuration interaction (CI) method, 
which is widely known in quantum chemistry and has recently been used to invest- 
igate finite graphene quantum dots^'*^. This multi-determinantal approach, built in 
molecular orbital space, provides excellent estimates of the ground and excited 
states^^. To avoid the exploding size of the fuU-CI space with increasing width, we map 
the whole system into a complete active space (CAS) with the same number of valence 
and conduction bands around the Fermi energy at each /c-point and consider aU 
possible electron occupancies to construct the many-body basis. In our calculations, 
we consider 4 valence and 4 conduction bands around the Fermi energy. This choice 
gives 4,900 many-body basis states in CAS for the half-filled case with S^"" — 0. Note 
that in the case of periodic systems, the scattering processes within different fc-points 
become crucial in determining the electronic or optical excitations. Although these 
scattering processes have been neglected in our study for simplicity, we confirm the 
qualitative observations by fmite size calculations that inherently take into account 
such scattering events. The hole doping has been performed by the systematic 
removal of electrons from CAS at each /c-point. This approach resembles the scat- 
tering process in angle -resolved photoemission spectroscopy (ARPES). With hole 
doping, the number of many-body basis states in CAS gradually decreases. Note that 
electron doping is expected to show simliar results, since the electron-hole symmetry 
is conserved in the nearest-neighbour Hubbard model. 



SCIENTIFIC REPORTS | 2 :519 | DOI: 10.1038/srep00519 



6 




//// 



Figure 7 | Schematic representation of the possible appUcation of 
ZGNRs. (a) The semiconducting and antiferromagnetic ground state of 
ZGNR with antiparallel spin orientation between opposite edges in 
absence of gate bias (V^). (b) Hole doping due to the application of gate 
bias makes the ZGNR magnetic and conducting serving as spin filter. 
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